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Abstract. 

We review the main results from recent numerical simulations of turbulent frag- 
mentation and star formation. Specifically, we discuss the observed scaling relation- 
ships, the "quiescent" (subsonic) nature of many star-forming cores, their energy 
, balance, their synthesized polarized dust emission, the ages of stars associated with 

' the molecular gas from which they have formed, the mass spectra of clumps, and 

the density and column density probability distribution function of the gas. We 
then give a critical discussion on recent attempts to explain and/or predict the star 
formation efficiency and the stellar initial mass function from the statistical nature 
of turbulent fields. Finally, it appears that turbulent fragmentation alone cannot 
^■f^ ■ account for the final stages of fragmentation: although the turbulent velocity field is 

able to produce filaments, the spatial distribution of cores in such filaments is better 
explained in terms of gravitational fragmentation. 
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1. Introduction 

It is well known that stars form in dense cores within molecular clouds, 
. but the physical processes that control the formation of low-mass stars 

i within these cores are not well understood yet. The "standard" scenario 

' assumes that low-mass protostellar cores are in quasi-static equilibrium 

supported against gravitational collapse by a combination of magnetic 
O I and thermal pressures (see, e.g., IShu et al.l HWH . In this scenario, 

• a core forms stars once magnetic support is lost through the pro- 

, cess known as "ambipolar diffusion": neutral gas particles slowly drift 

through the ions, which are held up by the magnetic field, allowing the 
core to eventually attain the "critical" mass-to-flux ratio, beyond which 
the gravitational energy exceeds the magnetic energy and collapse sets 
■ in from the inside-out. 

The more recent turbulent scenario, on the other hand, suggests 
that cores are formed by compressible motions in the turbulent ve- 
locity field of the cloud flSasao, 1973; 'Elmegreen, 1993| Padoan, 1995} 



Ballesteros-Paredes et al., 1999a^,Klessen et al., 2000( ). Those cores with 



an excess of gravitational energy collapse rapidly to form stars, while 
the others with sufficiently large internal or kinetic energies re-expand 

once the turbulent compression subsides ( [Taylor et al., 1996| Vazquez-Semadem et al., 2002 bl) 



Although the standard scenario has explained many observational 
features of protostellar cores, it does not address the problem of how the 
quasistatic cores can form and persist within a turbulent environment. 
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On the other hand, the turbulent scenario naturally incorporates the 
process of core formation, and is increasingly succesful in reproducing 
the observed core properties. In this paper, we will review the main 
results of the turbulent scenario, the observational facts that numerical 
works have successfully reproduced, the simulations' predictions con- 
cerning the star formation efficiency and the attempts to explain the 
initial mass function from statistical considerations of a turbulent field. 
Finally, we discuss the final stages of fragmentation, which appear to 
be gravitational rather than turbulent. 



2. Turbulent Fragmentation 



The process through which a chaotic velocity field produces a clumpy 
density structure in the gas within in a few dynamical timescales is 
referred to as "turbulent fragmentation". 

It is generally accepted that small-scale turbulence can provide sup- 
port against global gravitational collapse (e.g.. lChandrasekhar fc Fermi| 
I1953|> . In fact, if only small-scales modes are considered, and if they 
are more or less isotropic, they provide extra pressure. In this case, an 
"effective" Jeans length can be defined in terms of the "effective" sound 
speed, Cgg = + ^^^ms' where Urms is the turbulent velocity dispersion. 



which in principle depends on the scale ( Bonazzola et al., 1987 ). Thus, 
if the scale of a system is smaller than the modified Jeans length, it is 
said that turbulence inhibits collapse. This is valid only if turbulence 
is isotropic, and operates on small scales. In reality, interstellar tur- 
bulence is a multiscale phenomenon, and its largest scales are in fact 



comparable to the size of the system ( Ballesteros-Paredes et al., 1999at 



Brunt, 2003), and contain most of the kinetic energy. If this occurs, 
the large- and intermediate-scale modes of turbulence can produce 
fragmentation, while the small-scale modes can provide additional pres- 
sure. The structures produced by the fragmentation (cores) can either 
proceed to collapse, or redisperse and merge back into the ambient 
medium, depending on the balance between the energies within them 
( Sasao, 1973l|Elmegreen, 1993||Ballesteros-Paredes et al., 1999aD . Thus, 



supersonic turbulence plays a dual role, providing support for large 
structures while simultaneously inducing their fragmentation into smaller- 
scale subunits. In what follows, we review several predictions and com- 
parisons with observations from numerical simulations of turbulent 
fragmentation. 
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Turbulence decays fast. Molecular clouds (MC) are known to be tur- 



bulent, with motions that are supersonic at scales > 0.1 pc ( Larson, 1981 1 
Early concerns that turbulence should be dissipated rapidly in shocks 
( Goldreich &: Kwan, 1974 1 were shortly afterwards argued to not be 



applicable if turbulence consisted of MHD waves ( Arons &: Max, 1975 1 



However, an important result found in numerical simulations is that 
the original concerns were well founded, and that strong MHD tur- 
bulence decays rapidly even in the presence of strong magnetic fields 
(jPadoan fc Mordlun dj I Mac Low ct ai.; "TW; 'Ston e et al.l mM 

iMac Low. .1999:. A vila- Reese &: Vazquez-Semadeni^.2001. see, however, 
ICho et ahl I2UU2I for an alternative viewpoint). Thus, it appears that 
MC turbulence must be continually driven during the entire lifetime 
of the clouds (which, on the other hand, probably does not exceed 
several Myr - iBallestcros-Pa redes et al.l 1999b). The driving mecha- 
nism, however, is probably not restricted to the stellar activity internal 
to the clouds. Turbulence in clouds may be generated by the very 
compressive process that forms the cloud, through bending-mode insta- 
bilities aided by c ooling golini fc Walder, 1998; Klein fc Woods, 19981 



Koyama Hi Inutsuka, 2002 1, as proposed bv iBallesteros-ParedesTFaLl 



(|1999a|l and lVazquez-Semadeni, et al 



Energy balance. An evident result from numerical simulations is that 
the structures formed in those models evolve dynamically because they 
are part of a globally supersonic turbulent flow. Then, turbulence con- 
tinuously reshapes their structure, producing changes of mass, momen- 
tum and energy in them ( Ballesteros-Paredes et al., 1999a I . Clouds in 



simulations are found to be roughly in energy equipartition, but not in 
virial equilibrium, since the time derivatives of the moment of inertia 
and the surface terms in their virial balance are at least as important as 

the v olume energy terms^ ([B allesteros-Pared es Hi Vazquez-Semadeni} 

119951 [Ballesteros-Paredes fc Vazquez-Semadenilll997l seelBallesteros-Paredesl 



120031 for a review) . That the surface terms are large implies the ex- 
istence of strong fluxes of all physical quantities (mass, energy and 
momentum) across cloud "boundaries" which, on the other hand, are 
rather arbitrary. 



^ This does not means that "Virial Mass" estimates are wrong. In fact, those 
estimations are based on equipartition between gravitational and kinetic energies 
which, as mentioned above, are roughly in equipartition. Instead, "Virial Mass" 
estimations should be termed "Equipartition Mass" estimations. 
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Larson's Relations. Numerical simulations including turbulence have 
been helpful to understand the origin of the observed scaling relation- 
ships, Av oc i?^/^ and (p) oc ( [Larson, 1981 ). It is more or less 



accepted that the first relationship is a direct consequence of the tur- 
bulent power spectrum of a field dominated by shocks. However, it has 
been suggested that the second relationship may be an artifact of the 
observations ( Kegel, 1989| Scalo, 1990 1. In fact, it has been shown from 
numerical simulations that the relationship does not exist in real space 



( Vazquez-Semadeni et"aLlll997|IBanesteros-Paredes fc Mac Low|l2nn2|> . 

but it appears in the observational (position-position-velocity) space as 
a consequence of the thresholding done in the process of defining clouds 
or clumps (Ballesteros-Paredes Mac Low, 2002). 



Fast cloud and star formation. If, as suggested by simulations, clouds 
are far from virial equilibrium and supersonic turbulence decays fast, 
then probably cloud formation can be faster than earlier estimates, 
which were ~ 20-30 Myr (.Blitz fc Shu, 1980D. This ag rees w ith the 
suggestions by, among othersTH'asaol()1973|).IHunter et al.l()1986|),|Elmegreen 



()1993|) . IBallesteros-Paredes et al.1 (11999a ) , that clouds are turbulent 
density fluctuations within the global ISM flow. 

iBallestero s-P aredes et aP (1999b) have shown that the global tur- 
bulence driven by star formation events can produce locally gravita- 
tionally unstable clouds with characteristics similar to those found 
in Taurus. They argue that rapid assembly of molecular clouds by 
large-scale streams is a very plausible mechanism to form clouds, in 
which, furthermore, stars older than ~ 5 Myr are not found^. Indeed, 
IHartmann et al.l (|2()()H1 tabulate the ages of stars in 13 nearby star- 
forming regions. For those regions with stars older than ~ 5 Myr, there 
is no molecular gas associated, suggesting that the time scales for both 
cloud- and star- formation are short. 



2 Note that lPaila fc Stahierl f2000. l2002ll argue that nearby MC have been form- 
ing stars in the last 10 Myr or more, with a recent burst of star formation. However, 
IHartmannl J^HH^ argues that their conclusions are skewed by a statistically small 
sample of stars with masses larger than IMq, and by biases in their birthline age 
corrections. He notes that the picture bv lPalla fc Stahlerl 120021 fal requires the last 
1-2 Myr to be a special epoch for most MCs; (b) implies that most MCs are forming 
stars at extremely low rates, if any; and that (c) the apparently oldest stars are 
systematically higher in mass, implying that for most of a typical MCs star-forming 
history, the Initial Mass Function was strongly skewed. Item (a) seems implausible, 
and items (b) and (c) are contradicted by observations (see |Bartmann, 120031 and 
references therein). A more plausible explanation of the observations, he argues, is 
that the "tail" of older stars is really the result of including older foreground stars, 
as well as problems with the isochrone calibration in the higher mass stars. 
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Isolated and clustered modes of star formation. In numerical simu- 
lations, the structure of a turbulent molecular cloud depends on the 
largest scale of the turbulent motions ( Klessen et al., 2000} Heitsch et al., 2001a I 



In particular, for turbulence driven at large scales, "clustered" col- 
lapse occurs predominantly in filaments in a more or less efficient 
way, while for small-scale driven turbulence, a more "isolated" mode 
of star formation occurs (see IMac Low fc Klessen| 120031 for a review) . 
However, a potentionally confusing point is worth noting here. For the 
star formation community, a clustered region is one in which many 
stars are formed in very compact regions. For instance, while the Orion 
nebula has a star density of ~ 10^ — 10^ stars per cubic pc, the Taurus 
Molecular cloud contains only about 200 T Tauri stars. Thus, Orion is 
conventionally considered to belong to the clustered mode, and Taurus 
to the isolated one. However, from the point of view that Taurus is a 
coherent filamentary region, apparently due to an external compres- 



sion ( Ballesteros-Paredes et al., 1999b) that has formed stars almost 



simultaneously, it would belong to the "clustered" mode as defined by 
IKlessen et aD(|M7ni) . 

Existence of Quiescent Cores. Observational evidence suggests that 
low-mass stars form from molecular cloud cores whose column density 
profiles may resemble Bonnor-Ebert ( Ebert, 1955| Bonnor, 1956D equi- 



librium spheres, and whose velocity dispersions are transonic, or even 
subsonic JMyers, 1983l|Jijina et al., 1999t[Casem et al., 2002] ). For this 



reason such cores are usually termed "quiescent". In the standard 



scenario of star formation ( |Shu et al., 1987 1 quiescent structures are 



naturally explained as consequences of the quasistatic contraction pro- 
cess. In the turbulent scenario, however, protostellar cores are transient 
features generated by the dynamical fiow in the cloud. In a recent pa- 
per, iBallest eros-Paredes et al. (200_3 ) demonstrated that indeed, tran- 
sient, dynamic cores have an angle-averaged column density structure 
that often resembles hydrostatic Bonnor-Ebert profiles, such as ob- 
served cores (|Alves et al., 2001| [Evans et al., 2001] [Shirley et al., 2000} 



Johnstone et al., 2000D . They argue that this kind of analysis in obser- 



vational work is not an unambiguous test of hydrostatic equilibrium. 

Furthermore, IKlessen et al.l ()2003|1 investigate the velocity structure 
of these cores. They exhibit sub- and transonic velocity dispersions 
and, in self-gravitating cases are in approximate equipartition between 
the kinetic and gravitational energies. They thus closely resemble the 
observed population of so-called quiescent low- mass cores. The fact 
that dynamically evolving cores in highly supersonic turbulent fiows can 
have subsonic velocity dispersions is a direct consequence of a turbulent 
energy spectrum that decreases with decreasing length scale. Note that 
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this also implies that the quiescent nature of the cores is not necessar- 
ily indicative that they constitute the dissipative scale of turbulence, 
as is sometimes suggested (e.g., IGoodman d^rmi) . The observed 
properties of quiescent low-mass protostellar cores in molecular clouds 
thus do neither directly imply the existence of strong magnetic fields for 
support against collapse nor quasistatic evolution through ambipolar 
diffusion. 



Polarized thermal dust emission. Other studies have focused on the 
polarized emission of dust from numerical simulations. The main results 
are, first, that cores in the magnetized simulations exhibit degrees 
of polarization between 1 and 10%, regardless of whether the tur- 
bulence is sub- or super- Alfvenic ( Padoan et al., 2001a ); second, that 
submillimeter polarization maps of quiescent cores do not map the mag- 
netic fields inside the cores for visual extinctions larger than Ay ~ 3 
( Padoan et al., 2001a I; third, that the Chandrasekhar- Fermi method of 
estimating the mean magnetic field in a turbulent medium underesti- 
mates very weak magnetic fields by factors up to 150 C Heitsch et al., 2001bD ; 
and fourth, that limited telescope resolution leads to a systematic over- 
estimation of the field strengths, as well as to an impression of more 
regular, structureless magnetic fields (Heitsch et al., 2001b I. 



Density and column density PDF of turbulent gas. Several numeri- 
cal studies have investigated the statistical distribution of the density 
field. The volume- (mass-) weighted probability distribution function 
(PDF) of the density f{p)dp is defined as the fractional volume (mass) 
occupied by gas in the density range (p, p + dp) . 

It has been shown that, for a polytropic gas (i.e., for a gas where P oc 
p'^) the volume-weighted density-PDF depends on the effective poly- 
tropic exponent 7 of the fiuid. In particular, for an isothermal flow, (7 = 
1), the PDF is lognormal (IVazquez-Semadeni, 1994HPadoan et al., 1997| 
Passot fc Vazquez-SemaJeni, 1998||Scalo et al., 1998MOstriker et al., 2001} 
LieT^I', 2003 1, because the density fluctuations are the result of a 
stochastic multiplicative process, and so the distribution of its loga- 
rithm is Gaussian according to the Central Limit Theorem (| Passot fc Vazquez-Semadeni, 1998 ) . 
In the more general polytropic case (7 ^ 1), [Passot &: Vazquez^^emanem| 
H1998|) extended this result and showed that the PDF should approach 
a power law at high (resp. low) densities for 7 < 1 (resp. 7 > 1). 
Subsequently, other workers seem to have found somewhat different 
numerical results. FWada &: Norman1H2001() claim that a lognormal den- 
sity PDF also appears in their non-isothermal simulations of galactic 
disks, while ILi et al.1 (|2003|1 report that the shape of the PDF may 
depend on the scale in which the turbulence is driven (see also lKlessen| 
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I2UUUI However, the discrepancy may be apparent. The lognormal PDF 
of IWada &: Normaiil H2L)U1() appears only at the highest densities in 
their simulations, where they are nearly isothermal, while the result 
of ILi et al.l (|2()(),Sj) may be due to a different choice of the turbulent 
driving, as already noted by Passot k. Vazquez-Semadem (|1998j) and 
bvlKlessenll 



Finally, Passot &: Vazquez-Semadem (|2()(K-ij) have investigated the 
density PDF in the magnetic isothermal case, together with the role of 
the magnetic pressure, finding that the latter is only poorly correlated 
with the density, because of the different scalings of the magnetic field 
strength with the density fluctuations for the different types of nonlin- 
ear MHD waves. This implies that, rather than providing a restoring 
force, the magnetic pressure gradient acts more as a random forcing, 
and therefore in general has little effect on the density PDF, which in 
most cases retains its lognormal shape in the isothermal also 
found by lOstriker et aJ.1 pHTT) . 

Since direct measurements of the density field in the ISM are not 
available, it is important to understand the PDF of the column density 
field. For that purpose, Ostriker ct al. ( 2001 ) measured the PDF of the 
column density for numerical simulations, finding that this distribu- 
tion has a similar shape as that of the underlying density field, but 
with smaller mean and width. However, Vazquez-Semadeni &: Garcia 



H2001() argued that if the events are separated by distances larger than 
a correlation length, then the distribution of the density field should 
approach a Gaussian shape, although the convergence is very slow in 
the high density tail, so that the column density PDF is not expected 
to exhibit a unique functional shape, but to transit instead from a 
lognormal to a Gaussian form, as the ratio of the column length to the 
decorrelation length increases. This result, then, explains the nearly 
lognormal distribution found by lOstriker et al.l I 



IB lit z fc Williams (199^ performed a comparison between column 
density PDFs from models and observations by analyzing the PDF of 
the antenna temperature Ta of ^^CO for Taurus, assuming that Ta is 
proportional to the column density of the cloud. In their Fig. 1, it is 
clear that the central region behaves closely to a lognormal distribution. 
However, at larger Ta the distribution found for Taurus has an extended 
tail. This result may be due to the fact that in Taurus self-gravity has 
produced high-density collapsing cores. This is in fact apparent in the 
evolution of the volumetric density PDF shown by IKlessenI (|2000|l (see 
also Fig. 9 in ILi et al.l I2003|l , it is clear that self-gravity allows for the 
production of such tails in the volumetric density PDF. Thus, such tails 
should appear also in the column-density PDF. 
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Mass Spectra. The mass distribution of clouds and clumps in turbu- 
lent simulatio ns has been studied by several workers dVazquez-Semadeni et al., 1997) 
Klessen, 200T| [Ballesteros-Paredes Mac Low, 2002||Gammie et al., 2003D . 



Their results generally give a mass spectrum of the form dN/ dM oc M'^ 
with —2.5^a^ — 1.5. Since some of the results are discussed in the 
chapter by Klessen (this volume), here we just comment that several 
studies have recently tried to relate the mass spectrum of compact cores 
with the slope of the mass distribution of the newborn stars (IMF). 
However, this appears questionable because (a) the mass distribution 



of cores changes with time as clumps merge together ( Klessen, 2001 1; 
(b) cores generally produce a cluster of stars, not a single one, and 
so the relation between the masses of cores and of individual stars is 
unclear, and (c) many other processes are at play in determining a star's 
final mass, such as competitive accretion (see, e.g.. Bate, this volume) 
and stellar feedback through winds and outflows. In other words, a core 
with mass M will not necessarily form, in general, a star of the same 
mass, or even proportional to it. 

The Star Formation Efficiency. Studies bv lPadoanI p995jl and lKlessen et aP 
(f2000.} have proposed that turbulent fragmentation has an important 
role in determining the Star Formation Efficiency (SFE). Moreover, 
the latter authors found that the efficiency in numerical simulations 
(measured as the fraction of mass in collapsed objects) increases as the 
energy injection scale is increased. 

As remarked by IPadoanI (|1995j) . the compressible cascade regime 
must end at a sonic scale, so that the turbulence within structures 
smaller than this size is subsonic and thus cannot produce any further 
subfragmentation. Thus, in this scenario, subsonic cores are the natural 
end of the compressible part of the cascade in the ISM^. Moreover, 
when turbulence becomes subsonic, it ceases to be a stronger source 
of global support than thermal pressure, and the Jeans criterion can 
be applied to determine whether a core is gravitationally unstable 



and should proceed to collapse ( Padoan, 1995|), o r else rebound and 



merge back into its environment (^Imegreen, 1993' 'Taylor et al., 1996} 



Ballesteros-Paredes et al., 1999aj.Vazquez-Semadeni et al., 2002b |. This 



effectively sets an upper limit to the SFE because if the cloud is glob- 
ally supported by its internal turbulence, only massive enough clumps 
among those with sizes A smaller than the typical scale where the turbu- 
lence becomes subsonic, are susceptible of local collapse ([Padoan, 1995 1 . 



^ This does not necessarily imply that the dissipative (or "inner") scale of the 
turbulence has been reached, but only that at this point, the cascade can turn into 
an incompressible one. Similarly, it does not imply that there is a unique sonic scale. 
A large scatter in the velocity dispersion of structures with a given size is observed in 
simulations (e.g., |Vazquez-Semadeni et al.| I19t?71 [Vazquez-Semadeni, et ar]l2UU3ll . 
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In this sense, star-forming cores constitute the intersection of the set of 
subsonic regions with the high-mass, super-Jeans tail of the mass distri- 
bution ([P adoan, 1995 : Ballesteros-Paredes &: Vazquez-Semad eni, 1997} 
Padoan fc Nordlund, 2002,, , Vazquez-Semadeni, et al., 2003, ) . In partic- 
ular. IPadoanI (fl995j) gave a semiphenomenological theory which pre- 
dicted the SF efficiency and mass distribution of protostellar cores in 
turbulent MCs from the turbulent flow parameters, using the Larson 
scaling relations as empirical input. In that work, however, he used 
a questionable transformation between density and Jeans mass (see 
IScalo et al.l 11998 for a discussion). 

In this scenario, since progressively smaller density peaks contain 
progressively smaller fractions of the mass, one can expect the SFE to 
decrease with decreasing sonic scale, the latter being a function of the 



turbulent flow parameters (Vazquez-Semadeni, et al., 2003 >. Indeed 



the latter authors found that the star formation efficiency in numerical 
simulations scales with as e"''*''/'^'' , where Aq is a reference scale, 
which for the set of simulations fitted, has a value of ~ 0.02 pc. 

The Initial Mass Function from turbulent fragmentation. Interesting 
attempts to explain the Initial Mass Function (IMF) from turbulent 
models have been performed bv lPadoan et al.l(|1997|) and lPadoan fc Nordlundl 
(|2()()2j) . In particular, in the latter work, the authors propose a model in 
which the power spectrum of the turbulence is a power law, the cores 
are formed by shocks, and the typical size of a dense core scales as 
the thickness of the postshock gas; then the number of gravitationally 
unstable cores is given by 

N{m)dlogm (X m~^^^^~^^ p{mj)dm.j^dlogm, (1) 

where /? is the exponent of the energy power spectrum of the turbulence, 
E{k) oc k~^, and the integral of p{mj) is the fraction of cores of mass m 
with gravitational energy larger than the thermal energy. In particular, 
(Padoan &: Nordlund, 200*2 ) argue that the value of (3 = 1.8, which is 



close to the value found in their numerical simulations, and is close to 
the average value inferred from observations of the velocity dispersion- 
size relationship, gives a mass dependence in eq. (^ consistent with the 
Salpeter IMF. 

Caveats in the derivation of the SFE and the IMF. As interesting and 
promising as it is, the model by IPadoan &: Nordlundl (120021 is based 
on strong assumptions which we briefly discuss here in order to asses 
their feasibility. First, the number of cores as a function of scale, N{L) 
(their eq. [17]) is derived from scaling considerations as in Elmegreen] 



(|1997j) . Thus, it is independent of estimates of core properties (eqs. 
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[l]-[3], [11] in IPadoan fc Nordlundl IMl^ . while ideally, they should 
be self-consistent. Second, the derivation of estimate of N{L) is based 
on assuming that a large-scale (L2) simulation can be thought of as 
made of {L2/Lif space- filling simulations of scale Li, thus not taking 
into account intermittency of the density and velocity gradient fields. 
Instead, it assumes that the small scales are uniformly distributed in 
space without belonging to localized larger-scale ones. Thus, it will 
be probably neccessary to do a correction analogous to 
(|1962j) or /3-model intermittency corrections to original 



Kolmogorov 



Kolmogorov 



(|1941j) model. This might alter the scaling. Third, their eq. (24), which 
gives the number of collapsing cores, assumes that all collapsing cores 
are subsonic and indivisible by turbulence. Specifically, if a core with 
mass mi < 1112 is contained within a core of mass m2, their eq. (24) 
will predict two stars of masses oc mi and oc m2, implying that the 
same mass used to form the first star, is also used to form the second 
one. Moreover, recent numerical simulations of the collapse of molecular 
cloud cores by (Goodwin et al., 20031 have found that the number and 
properties of objects formed in an individual core depends in a chaotic 
way on the details of the initial velocity field (assumed subsonic, with 
Mach ~ 0.3), showing that even subsonic cores fragment into a random 
number of stars. Fourth, in their treatment, competitive accretion and 
stellar feedback is not considered as agents that might change in a 
substantial way the IMF. However, the interaction of protostellar cores 
and their competition for mass growth from their surroundings are 
important processes shaping the distribution of final star properties 
(see the reviews bv fMac Low fc Klessen ( 200,'T)). and lL^on (2003.)). 

Finally, it is worth mentioning that in their treatment, both lPadoan &: NordlunHI 
H2002() and [Vazquez-Semadeni, et al. H2003() consider that the collaps- 
ing regions are those with subsonic turbulence and masses larger than 
the Jeans mass. However, there are examples in our Galaxy of compact 
(i? < 0.1 pc), massive (M ~ IQ'^Mq) collapsing cores with velocity dis- 
persions of the order of 5v ~ 4-12 km/sec (see, e.g.. lKurtz et aHl2()()()l 
IMoliTiari et"ail ^MM IFa.undez etaTl VBW\f . Although it is likely that 
they might suffer from internal turbulent fragmentation, it is not clear 
whether there will be one or several (and how many) sites of collapse in 
them. In fact, ~ 20% of the objects observed bv lFaundez et al.1 ( 



(e.g., IRAS 16272-4837) seem to be good candidates for isolated, mas- 
sive star formation regions. Thus, these still need to be accommodated 
by the emerging theory of turbulent star-formation. 

Gravitational fragmentation on a filament produced by large-scale streams 
One important question is whether there is again a transition from 
turbulent fragmentation to gravitational fragmentation at some small 
scale, and if so, what is that scale. 
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In the particular case of Taurus, once turbulence has induced the for- 



mation of the filamentary structure of the gas ( |Ballesteros-Paredes et al., 1999bD , 
gravitational fragmentation in the filaments should be considered as a 
mechanism for forming cores and stars. In fact, following lLarsonl p985)l . 
ISartm ann (2002) has shown that the Jeans length and the free-fall 
timescale for a filament, are given by 



Xj = 1.5 Tio V PC, 
T ~ 3.7 T^/^ Ay^ Myr. 



(2.1) 
(2.2) 



where Tiq is the temperature in units of 10 K, and Ay is the visual 
extinction through the center of the filament. By using a mean visual 
extinction of starless cores Av- ~ 5 dUnishi et al., 1998D , eq. (EU gives 
a characteristic Jeans length of Aj ~ 0.3 pc, and collapse should occur 
in about 0.74 Mvr. iHartmannl (|2002|1 argues that, since the number of 
young stellar objects per length scale in Taurus is about 3 — 4 pc~^, 
then self-gravity is working and producing the cores and stars. 



^^^""^ T = 0.1 Myr 1=0.2 Myr 


t = 0.3 Myr 
-Aj 


L t = 0,4 Myr t = 0.5 Myr t = 0.6 Myr 
l^^j -Aj -Aj 


t = 0.7 Myr 



Figure 1. Evolution of the column density of an SPH simulation (courtesy of R. 
Klessen). The filament in the first frame (before self-gravity is turned-on) shows that 
turbulence is responsible for forming this kind of structures. The small bar in the 
bottom-left of each frame denotes the Jeans length (eq. 12.11 ) at this time. At later 
times, self-gravity is turned on and the filament suffers gravitational fragmentation 
on a free-fall timescale (eq. 12.21 '). The size of each box is 6 pc. 

In order to test this idea in the simulations, we analyze an SPH 
calculation performed by R. Klessen (see his contribution in this volume 
for more details). This simulation has been performed without gravity 
until a particular, well defined elongated structure is formed. At this 
time, self-gravity is turned on. The results are in good agreement with 
IHartmannl H2002|) . In Fig. ^ we show eight column density frames of 
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the simulation. The first frame shows the structure just before self- 
gravity is turned-on. We see that the filament forms cores in a fraction 
of a Myr. The timestep between frames is 0.1 Myr. The mean surface 
density for the filament is 0.033 gr cm~^, corresponding to a visual 
extinction of ~ 7.5. Using eqs. (|2.1j) and (|2.2|) . this value gives a Jeans 
length ofAj~0.2pc, and a collapsing timescale of r ~ 0.5 Myr. Note 
from Fig.^that the first cores appear roughly at r ~ 0.3 Myr, although 
the final structure of collapsed objects is clearly defined at t = 0.5 Myr. 
Moreover, the typical separation between the cores corresponds roughly 
to the Jeans length, which is indicated by a small bar in the figures. 

We then conclude that even though the filament is not in hydrostatic 
equilibrium, gravitational fragmentation may produce the structures 
seen either in the simulations, and in actual clouds. This means that 
turbulence causes first the production of the filament, and then self- 
gravity takes over to form the collapsed objects. However, since (a) 
the filament is not a perfect cylinder, (b) the collapsed objects are not 
perfectly equally spaced, and (c) there is a range in timesteps in which 
the cores are formed (between t ~ 0.3 and 0.6 Myr), it becomes clear 
that the final properties of the star forming regions still depend on the 
initial conditions induced by the global turbulence. In a similar way, the 
final mass of the stars may still depend on the accretion history from 
the surrounding cloud material (see chapter by Klessen, this volume). 

4. Conclusions 

This review shows that numerical simulations of turbulent fragmen- 
tation have been quite successful in reproducing a number of observa- 
tional properties of molecular cloud cores, while simultaneously making 
a number of predictions about the behavior of turbulence in molecular 
clouds and their cores, as well as providing an explanation for another 
set of observed features. Although much work remains to be done, and a 
full theory is still elusive, the progress made is highly encouraging, and 
makes the dynamical picture of star formation a plausible alternative to 
the ambipolar-diffusion-mediated magnetostatic-equilibrium standard 
model. 
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